im = imread('kofka_ring.tif');
[Ix,Iy] = ImageDerivatives(double(im));
lap = abs(Deriv2Laplace(Ix,Iy));

for i = 1:2:20
    i
    temp = lap;
    temp(temp<i) = 0;
    imshow(temp);
    pause
end